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d' 

O^ ' We make a comparative study of quadrature squeezing, photon-number distribution and 

^ ' Wigner function in a decayed quantum system. Specifically, for a field mode prepared initially 

in cat states interacting with a zero-temperature environment, we show that the rate of 
reduction of the nonclassical effects in this system is proportional to the occurrence of the 
decoherence process. 

PACS numbers: 42.50.Dv,42.50.-p 

^, 

qh; I. INTRODUCTION 

^ . Decoherence process represents the transformation of superposition states into statistical mix- 

t^ ' 

Q \ ture states, i.e. the off-diagonal elements of the system are suppressed. This can occur through, 

\0 '• 

^f^ • e.g., the interaction between system and environment. Actually, the decoherence process is impor- 



o\ 



r^ I tant not only for understanding the quantum-classical transition 



but also it may eventually be 
useful for applications that require keeping coherence in mesoscopic or macroscopic systems, such 



as quantum com 



lutation (^. Furthermore, the decoherence is at the heart of the quantum theory 



^ ■ of measurement 

H : 

5t , On the other hand, superposition principle is at the heart of the quantum mechanics. It implies 

that probability densities of observable quantities usually exhibit interference effects instead of 
simply being added. The most significant examples reflecting the power of such principle are 
the Schrodinger cat states [J], which exhibit various nonclassical effects, such as squeezing, sub- 
Poissonian statistics and oscillations in photon- number distribution 5|-l7|, even if their components 
are close to the classical ones [8i],i.e., they are minimum uncertainty states and exhibit Poissonian 
distribution. These states can be defined as: 



a 



,^ = A2[|a)+exp(z0)|-a)], (1) 



where |a) is a coherent state with complex amphtude a , (/> is a relative phase and A is the 
normalization constant having the form 

^ = (2) 

2[l + exp(-2|a|2)cos(/)]' ^ ^ 

Specifically, for cj) = 0,-k and tt/2 state ([1]) reduces to even coherent (ECS), odd coherent (OCS) 
and Yurke-Stoler (YSS) states, respectively. It is worth mentioning that there are two regimes 
controlling the behavior of the states ( [1]), which are microscopic regime for small values of |a| 
(i.e., when the "distance" between the components of the cat is small) and macroscopic regime 
for large values of |a| [9]. In fact, these states are more nonclassical in the microscopic regime. 
In other words, the amount of nonclassical effects, such as the negative values in Wigner function 
and the oscillatory behavior in the photon-number distribution for these types of states are more 
pronounced in the microscopic regime than in the macroscopic regime. For more details about 
states ([1]), such as their generations and their properties when they are evolving in various optical 
systems, one can consult the review article 10|, and references therein. 

In this article we study the relation between the decoherence process and the occurrence of the 
nonclassical effects in a decayed quantum system. More precisely, we compare development of non- 
classical effects in both quadrature squeezing and photon-number distribution with the occurrence 
of interference pattern in the Wigner function. We perform such a comparison for the field mode 
prepared initially in the state ([T]) (described by the density matrix p) which interacts with zero- 
temperature environment. The master equation in the Born-Markov approximation describing the 
system is |ll| 



-Q- = -{2apa^ -a^ap- pa^a), (3) 

where 7 is the decay constant and a [a') is the annihilation (creation) operator designated to the 



mode of the field. The well-known time dependent solution for ([3]) is 12] 



P{t) = ^ X] exp(i(/>j.,/)(aj|aj/)^ ^'\y/Jlaj){^af\, (4) 

where ai = a, 02 = —a , p = exp(— r), r = ^7 is the scaled decaying parameter and 
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FIG. 1: Squeezing factor S2{t) of ECS case against the scaled decaying parameter r for a = 2 (solid curve), 
1.5 (star-centered curve), 1 (short-dashed curve) and 0.5 (long-dashed curve). 

II. QUADRATURE SQUEEZING 

As is vi^ell known squeezing is one of the most important phenomena in quantum optics because 
of its applications in various areas, e.g., in optics communication, quantum information theory, etc. 
[l3|. Squeezed light can be measured by a homodyne detection where the signal is superimposed 
on a strong coherent beam of the local oscillator. 

Here we investigate quadrature squeezing for the density matrix (jH). For this purpose we define 
the position and momentum operators, which are related to the conjugate electric and magnetic 
field operators E and H of electromagnetic waves, as 



X = -{a + a 



Y 



1 
2.' 



a — a') 



(6) 



where [X,Y] = §; then the uncertainty relation reads {{AX)^){{AYf) > -^ where {{AXf) = 
{X"^) — {X)"^ . Therefore, we can say that the mode is squeezed if Si{t) = 4((AX)^) — 1 < or 
S2{t) = 4((Ay)2) - 1< 0. 

Now squeezing factors Sj{t) for the system under consideration — restricting ourselves to ECS 
case — take the forms 



5i(t) = /i5i(0) 



4/xa^ 



-, S2{t) = ^iS2{Q) 



-Afia'^ exp(— 2a;^ 



(7) 



l + exp(-2a2)' '^ ' ^ ^^ ^ l + exp(-2a2) ' 

where Sj{0) are the initial squeezing factors. We have considered here a to be real. From these 



expressions it is clear that the quantum fluctuation of the field decreases exponentially as a result 
of its interaction with environment. More precisely, Sj{t) decay at the same rate as the intensity 



of the field [iJ]. Further, we see that squeezing exists provided that a and r are finite. Of course, 
the origin of these nonclassical effects is in the interference between the components of the cat. 
One can also check that when r is large enough, squeezing factors Sj{t),j = 1,2 tend to zero. 
In other words, the system tends to a steady state, which, in this case, is a pure state (vacuum 
state). In Fig. 1 we plot S2{t) for shown values of the parameters. Generally we can see that 
squeezing is more pronounced in the microscopic regime (when a is small). Furthermore, for a > 2 
the squeezing factor tends to zero regardless of the values of r. In this case the system becomes a 
statistical mixture state or a vacuum state if the values of r are small or large. This point will be 
clear when investigating the behavior of the Wigner function in section IV. 

III. PHOTON-NUMBER DISTRIBUTION 

Photon-number distribution P{n) is an integral part of the modern description of light, which 
can be measured by photon detectors based on the photoelectric effect. Further, one of the most 
interesting nonclassical effects emerging from the superposition principle is the oscillatory behavior 
in P{n). In general, such behavior is closely related to that of the Wigner function, however, this 
is a necessary but not sufficient condition. For example, the P{n) of ECS, OCS and YSS are 
completely different; whereas those of ECS and OCS exhibit pairwise oscillations in phase space 
(even number of photons can be observed for ECS and odd numbers for OCS), the distribution of 
YSS is a Poissonian even though the behavior of the Wigner function for these states is qualitatively 
similar. The second issue we want to address here is that in general the occurrence of squeezing 
in the quadrature variances does not need to be accompanied by oscillations in the P[n) and 
vice versa. For instance, binomial states [13] can exhibit quadrature squeezing even though their 
P{n) are close to Poissonian ones. In the same spirit, for ECS the oscillations in P{n) are more 
pronounced when the "distance" between the basis of the cat increases, however, this is not the 
case for the quadrature squeezing, which is completely suppressed for a > 2 (see Fig. 1). Now we 
investigate the sensitivity of the P{n) of the system under consideration to lossy mechanism. This 
quantity can be calculated easily {P{n) = {n\p{t)\n)) and one obtains 

P(n) = 2^ ^^"^ " exp(-^Q2) {1 + (-l)'^/(a) cos </>} , (8) 
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FIG. 2: P{n) of ECS case against n for a = 1 (a), 2 (b) and for r = (solid curve), 0.1 (short-dashed 
curve), 0.3 (long-dashed curve) and 1 (circle-centered curve). 



where 



f{a)=exp[-2a\l-f,)]. 



(9) 



By comparing the expression ([7]) with ([8]) we find that the dissipation is involved in two quantities by 
different ways and consequently the sensitivity of these quantities to lossy mechanism is completely 
different. As before, the origin of the nonclassical oscillations in the P(n) lies in the interference 
in phase space. Further, in ([8]) the interference term is decaying by the factor f{a) Hi and thus 
its contribution is more pronounced — oscillatory behavior can occur in P{n) — when a and /x are 
small. This situation is similar to that of the quadrature squeezing. We will discuss this point 
quantitatively in section IV by investigating the behavior of the factor /(a). In Figs. 2 we plot 
the P{n) for ECS against n for microscopic (a) and macroscopic (6) regimes, respectively, for given 
values of the parameters. By comparing the curves in Fig. 2a with those having the same values of 
r in Fig. 2b, one can conclude that the oscillations in P{n) for macroscopic regime are suppressed 
faster than those for microscopic regime provided that r is small. Also the comparison between 
the behavior of both the short-dashed curves in Fig. 1 and in Fig. 2a shows that the P(n) is more 
sensitive to dissipation than the quadrature squeezing is. This is clear as one can observe that the 
oscillations in the P{n) are completely suppressed, however, squeezing is still remarkable in the 
quadrature squeezing. The final remark is that the behavior of the P{n) when r = 0.3, 1 in Figs. 



2 is close to that for the statistical mixture of coherent states under the influence of the decay 
mechanism. 

IV. WIGNER FUNCTION 

Wigner (W) function is one of the quasiprobability functions, which carries full information 
about the quantum system. This function is sensitive to the interference in phase space and 
can be realized in optical homodyne tomography [la]. Here we use this function to study the 
decoherence of the system under discussion. The definition of the decoherence has been given in 
the Introduction. 

The W function can be defined as 

ww,t) = ^ J d^CewWC - l3*C)C''^\C,t), (10) 

where C*^"'^((",t) is the symmetrically ordered characteristic function having the form 

C(-)(C,t) =Tr[^(t)exp(atC-ar)], (11) 

where p{t) is the density matrix of the system, which for the system under consideration is given by 



B 



([!]). For the future purpose, we derive the W function following the same steps as in [8|. Thus we 
rewrite the W function in terms of the normally ordered moments of the creation and annihilation 
operators using the Baker-Hausdorff theorem. Therefore (jlOp takes the form 

w(M = ^t^-^^=W^i..., (12) 

TT^ ^-^ n\m\ 

n,m=0 

where we have used the abbreviation 

imn = /d2cexp(-i|cp + cp - cnn-cr 

(13) 

= (-i)"+"^a|^ / d\eM-m' + cp - en- 

Carrying out the integration in (J13p we obtain 

am+n 

I^^ = M-ir^-'g^,,,^ exp(-2|/?p). (14) 

After minor algebra and using the Rodrigues' formula for Laguerre polynomial, p4p reads 

/^n = 2"+V(-l)™r^"-'"^m!L:;-"(2|/3|2)exp(-2|/3|2), (15) 
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FIG. 3: The W function of ECS case for {a^r) = (1,0.3) (a); (a, r) = (2,0.3) (b). 

where L^{.) are the associated Laguerre polynomials of order m. 

On the other hand, the normally ordered expectation values {a)"^ (t)a"' [t]) associated with the 
density matrix (jl|) are given as |l7l |: 

2 

(a'f'^(t)a"(t)) =yl ^ e-x.v{i<^jj>){aj\aj,)a']a*^ ii"^ . (16) 

On substituting ([15]) and (fT6|) into (fT2|) we arrive at 



T^(/3,t) 



^exp(-2|/3|2) 



TT 



n,m=0 



^ exp(i(;^jj/)(aj|aj/)a"'a*r 



(1 -I ii.-i-//i 
^> \ ^, L;;-"^(2|/3p). 

(17) 
On using the generating function for Laguerre polynomials and the Taylor's expansion for the 
exponential function, (|17p reduces to the following closed form 



2A r 
VF(/3,t) = — exp(-2|/3-a^p) + exp(-2|/3 + a^|2) + 2/(a)cos(/)cos(4aV;uIm/3)exp(-2|/3p) 
vr L 

(18) 
where /(a) is given by ([9]). 

In general, the W function of ECS, OCS and YSS (at t = 0) are consisting of two Gaussian 
bells corresponding to statistical mixture of individual composite states and interference fringes in 
between originating from the superposition between different components of the states. Actually, 



these fringes represent the signature of the nonclassical effects. For this reason several articles have 
been devoted to deal with these fringes making them less or more pronounced by allowing the cat 



id ], and references therein). For 



states to evolve in different quantum optical systems (e.g., see 
the system under consideration we can easily conclude from (jlSp that as the interaction of the 
system with the environment is going on, the two Gaussian peaks of the statistical-mixture part 
move towards the origin and eventually emerge into each other. This is quite obvious since the 
centers of the peaks are exponentially decaying function of time. Furthermore, the amplitude of the 
oscillatory term goes down by the factor f{a) similar to the P{n). Such behavior can be explained 
as the flux of coherent energy transfers to the environment from the field and noise transfers to the 
field from the environment. More information about the system can be observed in Figs. 3a and b 
where we plot W{13 = x + iy) function for microscopic (a = 1) and macroscopic (a = 2) regimes, 
respectively. In both cases the scaled decaying parameter r = 0.3. From Fig. 3b it is clear that the 
optical cavity field tends to an approximate statistical mixture state, i.e., to a two-peak structure 
with negligible interference part. Actually the suppression of the nonclassical interference pattern 
in the W function does not mean that the system reaches its equilibrium states [1&\- Further for 
large interaction times the cavity field collapses to vacuum state irrespective of the type of the 
initial cat state. This can be checked from (jlSp as well as can be clearly seen in Figs. 8 and 9 in 
8|] (see curve 5 in these figures). This means that the superpositions of macroscopical cat states 
can be realized, but to have them surviving for some time the system must be completely isolated. 
Even a very slight interaction with the environment will very rapidly reduce the superpositions to 
the corresponding statistical mixture states. 

Now we turn our attention to the microscopic case (Fig. 3a). From this figure one can observe 
that the noise ellipse related to squeezed states is similar to that of squeezed vacuum states. 
The origin of this behavior is in the competition between the diagonal and off-diagonal elements 
of the system. Actually, in the microscopic regime the contributions of the statistical mixture 
components are located close to the origin of the phase space. Furthermore, the comparison 
between the behavior of quadrature squeezing, photon- number distribution and W function (i.e., 
the comparison of Fig. 1, Fig. 2 and Figs. 3a-b for the specified values of the parameters) shows 
that the occurrence of the nonclassical effects and decoherence phenomenon are qualitatively on 
the same level. More precisely, the more the system decoheres, the more the nonclassical effects 
decrease. This conclusion is completely different from that in [8]. The reason of this difference is 
that when the authors of [8j] compare the decay of the interference part in phase space based on 
W function (Fig. 6) with the behavior of the quadrature squeezing (Figs. 8, 9) they chose for the 
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FIG. 4: The function /(a) aganist a for t = 0.1 (solid curve), 0.3 (short-dashed curve), 0.8 (long-dashed 
curve) and 1.2 (star-centered curve). 

field amplitude a = 2, for which squeezing does not exist. Therefore both developments cannot 
be compared and thus they arrived at misleading conclusions. Furthermore, they explained their 
results by using series form for the W function (see ()17p : further in Eq. (5.13) of 8| there is a 
misprint in this expression, where no square root should be in its denominator) and concluded that 
" it is clearly seen that the Wigner function always decays faster than the second-order squeezing" . 
Actually, this discussion is not persuading because the expansion contains the terms of both the 
mixture and the interference components symmetrically. 

We conclude giving a quantitative analysis of the factor /(a) in Fig. 4. Such analysis can 
give insight into the occurrence (or nonoccurrence) of the decoherence process regarding to the 
values of a and the interaction time. As is clear from ([9]), /(a) exponentially decays whenever 
a increases provided that /i 7^ 1 (i.e., r ^ 0) and has its maximum value at // = (i.e., r is 
very large). Actually, Fig. 4, even if it is relatively simple, it can give the smallest values of a 
for which the system can be completely decohered for certain values of the interaction time. For 
instance, for r = 0.1, 0.3, 0.8 the corresponding smallest values are a = 5, 3, 2 for which the system 
is completely decohered. In this case the density matrix describing the system has typically the 
form Prit) = ^[\at){at\ + \ — at){—at\] where at = a^/Jl [7]. It is clear that these results agree 
with the fact that the nonclassical effects occur in the microscopic regime. Finally, it is worth 
mentioning that the decoherence in the present system can be overcome by including amplifying 
media in the cavity [l9|. 
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In conclusion, we have shown that the sensitivity of quadrature squeezing and W function to 
lossy mechanism is on the same level. This is not a surprising result since the W function is built 
on the complementarity of the canonical operators [20]. On the other hand, the P{n) is more 
sensitive to dissipation than the quadrature squeezing. Furthermore, the decoherence process is 
more visible in the macroscopic regime. Thus in a more realistic situation the generation and 
detection of a macroscopic superposition states is very difficult, due to the unavoidable coupling 
with environment and the consequent dissipation [21]]. Finally in the view of the quantities studied 
here, the nonclassical superposition states cannot be saved from decoherence. 
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